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Abstract. The Baryonic Tully-Fisher relation (BTF) can be substantially improved when considering that the galactic baryonic 
mass is likely to consist not only of the detected baryons, stars and gas, but also of a dark baryonic component proportional to 
the HI gas. The BTF relation is optimally improved when the HI mass is multiplied by a factor of about 3, but larger factors 
up to 1 1 - 16 still improve the fit over the original one using only the detected baryons. The strength of this improved relation 
is quantified with up-to-date statistical tests such as the Akaike Information Criterion or the Bayesian Information Criterion. In 
particular they allow us to show that supposing a variable M+/L ratio instead is much less significant. This result reinforces the 
suggestion made in several recent works that mass within galactic disks must be a multiple of the HI mass, and that galactic 
disks are substantially, but not necessarily fully, self-gravitating. 
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1. Introduction 



iMcOaugh et all fcOOlh (hereafter MSBD) have extended the 
Tully-Fisher relation llTullv & Fisher 1977) (TF) over 5 dex in 
stellar mass and 1 dex in velocity by correlating not the lumi- 
nosity with the rotation velocity, but the detected baryon mass 
with the rotation velocity. This is an important step toward un- 
derstanding the TF relation, because the physical causality link 
between mass and rotational velocity, gravitational dynamics, 
is much more direct than between stellar light production and 
rotational velocity. The traditional TF relation requires us to 
understand how the nuclear energy production inside stars may 
well be tightly correlated with the global rotation speed of the 
galaxy, while the baryonic TF (BTF) just requires finding the 
link between the baryonic mass and kinematics. 

The baryons considered by MSBD include the stars, the to- 
tal mass of which is inferred from a plausible stellar mass to 
light ratio, and neutral hydrogen HI from 21 cm measurements, 
including a cosmic helium fraction. Doing so, MSBD show 
with a sample of 243 galaxies that a scatter plot of velocity- 
mass looks straighter when summing the HI mass to the stel- 
lar mass. The argument advanced by MSBD is that a better 
relation exists between these baryons and rotational velocity 
than between light and rotational velocity because the low sur- 
face brightness galaxies in the sample have a very low stellar 
mass content. This makes sense because in the limit of a pure 
gaseous protogalactic disk no stars shine, yet the disk does ro- 
tate. 

Several reasons motivated us to reexamine this work. By 
naming the improved relation "baryonic" MSBD implies that 
the baryons are well represented by the detected stars and HI, 



while in many spirals often as much gas as the detected HI 
but in the form of H2 can be inferred from CO observations 
with the conventional proportiona lity X fac tor between the CO 
emission and the H2 mass ( lCombesl l2004). This point is how- 
ever not necessarily crucial because the CO emission is often 
roughly located in the optical disk, so the CO related H2 mass 
may be well absorbed, to a first approximation, in the mass to 
light ratio of the stellar mass. But for years now the universal- 
ity of the X fact or has been c hallenged. Recent rediscussions of 
this factor (e.g., BoulangeJ l2004l) favor an increased X in low 
metallicity, low excitation or cold conditions, typical of outer 
galactic HI disks, by at least one order of magnitude. 



Long ago Bosma (1981a) pointed out that the dynam- 
ical mass in spirals is known to be well correlated to a 
multiple of the detected HI (see also IHoekstra et aDl200ll) . 
So a natural hypothesis reinf orced by a series of argu- 
ments presented in Pfenniger, Combes & Martinetl i ll 994b . and 
Pfenniger & Combes ( 1994) as well as var i ants discussed 
by other authors (Henriksen & Widrow 1995; de Paolis et al. 
1995; Gerhar d & Silld ll996) is to assume that a sufficiently 
cold form of H2 can well make up a sizable fraction of the 
unseen baryonic mass in the disk. Moreover, several new theo- 
retical and observational arguments indicate t hat galacti c disks 
may be more massive than usually thought. Kalberla (2003) 
found that the mass distribution and the dynamics of the Milky 
Way is dominated by a dark matter disk. Massive disks are 
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also needed to explai n spiral structures in low surfac e bright 
ness galaxies (Fuchs 2003; Mayer & Wadslev 2004), as wel 
extended HI di sks like in the b 
]291 5 JMasset & Bureau! 2003) 
ties, Revaz & Pfenniger (2004) 



well 

as in exten ded HI disks like in the blue compact dwarf galaxy 
NGC291 5 JMasset & Bureaul|2003h . In the context of warped 
galaxies, Revaz & Pfenniger ( 2004J) have recently shown that 
the high number of warped spirals results naturally from bend- 
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ing instabilities, if the disk contributes to at least 60% of the 
total galactic mass, within 30 kpc. 

The true number of baryons contained in spiral galaxies can 
still be a multiple of the detected ones. Presently the hypothesis 
that MACHOS can represent an important part of the unseen 
baryonic mass in an extended halo seems to be disproved by the 
micro-lensing expe riments to wards the Magellanic clouds and 
the Galactic center dAlcock et al.l200lllAfonso et al.l2003t) . but 
a dark baryonic component closer to the outer disk, i.e., more 
closely related to the HI distribution, has not been investigated 
to a similar extent yet. 

If some sizable fraction of the mass correlates well with 
velocity, and the detected baryons within the radius of mea- 
sured velocities do represent a non-negligible amount of the 
total mass within this radius, it seems reasonable to check 
whether the full dynamical mass could not improve the cor- 
relation. Indeed, gravitational dynamics is little dependent on 
the nature of matter over a few dynamical time-scales. 

In order to investigate better the underlying apparent corre- 
lation, the analysis of the MSBD sample can easily be extended 
by up-to-date statistical tools. Below, we use statistical tools, 
mainly the Akaike Information Criterion (Akaike Il974h . and 
the Bayesian Information Criterion ( Schwarz 1978) that are not 
well known in astronomy, but have been often used in many 
other scientific fields over the last two decades. They provide 
interesting decision criterions for discriminating quantitatively 
between models using a different number of fitting parameters 
for the same data set. The curious reader will find in Liddle 
(2004) a short introduction for astronomers to these methods, 
and an application to the cosmological parameters. 

The purpose of the following Sections is therefore to extend 
the discussion about the BTF with more advanced statistical 
techniques than the ones currently popular in astronomy. 

2. Sample 

MSBD kindly provided us the sample that they used in their 
work, and to make a reasonable comparison we use here ex- 
actly the same sample, including the same assumptions about 
the mass-to-light ratios of the stellar and HI components 1 . The 
sample consists of 243 approximately face-on spirals with peak 
rotational velocities in the range 33 to 323kms~ 1 , almost a 
dex. For a more detailed description of the sample we refer 
the reader to the MSBD paper. 

3. Optimized HI scaling 

3.1. Full sample 

The first question is how scaling HI with a constant factor c dif- 
ferent from 1 changes the fit rms. Using a linear least-squares 
fit, for c in the range [0 - 15], we calculate the corresponding 
coefficients a and b in the least-squares relation linking the stel- 

1 In the following, and as commonly assumed in HI related works, 
by HI mass one means the 21 cm optically thin measured HI mass 
augmented by the cosmic He abundance. 
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Fig. 1. Rms of the linear least-squares fit as a function of the 
parameter c (see Eq.[0. The dotted line indicates MSBD's rms 
value (c = 1). 



lar mass M* inferred from the light and a constant M+/L ratio, 
the HI mass Mm, and the rotational velocity V c : 

log(M* + cM m ) « a + Z> log V c , ( 1 ) 

as well as the rms of the respective fits as a function of c, which 
are plotted in Fig. [2 Clearly, passing from c = (rms = 0.521) 
in the usual TF relation to c — 1 (rms = 0.329) in the BTF rela- 
tion decreases the fit rms by a factor of 1 .6, which is the main 
new point in MSBD's paper. Yet the curve rms(c) continues to 
decrease for c > 1 down to a shallow minimum at c = 2.98 
(rms = 0.3164). Beyond this optimal c, the rms grows slowly, 
such that up to c — 11.5 the rms fit is still better than for c = 1 . 
Fig. |2] shows the least-squares fit of the extended BTF (here- 
after EBTF) relation at the optimal value c opt = 2.98, where the 
linear fit parameters are a = 3.11 and b = 3.36. 

3.2. Comparing models 

With an additional parameter, we naturally expect that the rms 
of the EBTF model (Eq. [0 decreases with respect to the sim- 
pler BTF model (null hypothesis), where c remains fixed to 1 . 
Thus, even if the rms is lower, it does not mean that the data 
are better fitted by the EBTF model. 

To test this, we have used three different statistical tools. 
The first tool is classical, while the two more recent ones are 
less well known in astronomy, but have become popular else- 
where. The methods work under the assumption that the error 
distribution is Gaussian. A study of the rms distribution shows 
that the data points do verify reasonably well this assumption. 
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log V c 

Fig. 2. Least-squares fit of the Baryonic Tully-Fisher relation 
using Eq.Q]witha = 3.11, b = 3.36 and c = 2.98. 

3.2.1. Extra sum-of-squares F-Test 

This standard test consists of comparing the difference between 
the x 2 sums of two models, say model 1 and model 2, where 
model 1 is a simpler case of model 2 2 . Because model 1 has 
a lower number of parameters, we usually expect that x 2 \ is 
greater than x 2 2 with a difference corresponding to the differ- 
ence r of parameters between the two models. Thus, the ra- 
tio (x2\ - x 2 i) l r provides an estimate of the goodness of one 
model compared to the other. An F-value is defined by divid- 
ing the previous ratio by the normalised fraction x 2 2 l(N - ki), 
where N - k 2 is the degrees of freedom of model 2 with hi, its 
number of parameters. Explicitly, F is written as: 



F = 



(x 2 i ~X 2 i) I X 2 2 



(N - k 2 ) 



(2) 



Since F is a ratio of two x 2 distributions, it follows an F- 
distribution. We can then compute the probability P to find a 
value greater than F for such a distribution. If this probabil- 
ity is weak (typically less than 5%), we can conclude that the 
F-value obtained has a weak probability to result from the er- 
ror fluctuations and that the data is significantly better fitted by 
model 2. 

Assuming that the intrinsic errors cr,- are constant for each 
point i (cr, = cr), the x 2 terms in Eq. (|2j can be replaced by the 
sum-of-squares S 2(. 



S2 i = Y j (y j -f i (x j ))2, 



(3) 



2 By simpler case, we mean that model 1 is identical to model 2, 
with a number r of parameters fixed. The two models are also said to 
be nested. Since in MSBD's BTF model the parameter c is fixed to 1, 
this model is a simpler model than the EBTF one. 



where the yj correspond to the left hand side of Eq. Q and 
the fi(xj) to the right hand one. A rigorous desc ription of this 
classical test can be found, e.g., in lLuptonl l ll995t p. 98-101). 

Comparing the EBTF model with the simpler BTF model, 
with N = 243, k\ - 2 and r — 1, we obtain F = 19.9, to which 
the associated probability P is lower than 0.01%, which is the 
probability for the BTF model to be better than the EBTF one. 



3.2.2. Akaike's Information Criterion corrected (AlCc) 

The second performed test is derived form the Information 
theory, namely the Akaike's Informati on Criteri o n cor- 
rected (AICc) JSakamoto. Ishiguro & Kitagawal Il986t 
iBurnham & Anderson! I2004IK Since the formal derivation 
of this test is somewhat lengthy we refer the reader to the 
literature for more details (se e, e.g.. lBrockwell & Davisl F2002: 
IBurnham & Anderson! 120041) . However, its application is as 
simple as the Extra sum-of-squares F-Test, without the need 
to introduce additional assumptions 3 . The AICc consists of 
computing for the two models the values: 



2k i (k i + 1) 



(4) 



N-k-l 

where k, is the number of parameters of model i and 52; is the 
corresponding sum-of-squares defined by Eq. Q. In this for- 
mula, the first term is the entropy term, which usually decreases 
when the number of parameters ki increases. It is balanced by 
the second term linear in k t . The third term is the correction 
term especially relevant when N is small, so unimportant here. 
As for entropy in statistical mechanics, the value of depends 
on the choice of data unit, therefore its value has no absolute 
meaning, only the difference between two models has one: the 
relative information content. The best model is simply the one 
with the lowest J\ c value, where the decrease of the entropy 
term dominates. The difference Aj?l c = J\ c 2 - allows us to 
quantify how much a model is better than another. The proba- 
bility P24 that model 2 is better than model 1 is given by: 

= —■ (5) 



1 +e- 



The comparison of the BTF model with the EBTF one gives 
respectively m\ = -533.79 and J[ c 2 = -551.14. The so called 
"information ratio" fi.ifPi,! tells us that the EBTF model is 
5854 times more likely to be correct than the BTF one. 



3.2.3. Bayesian Information Criterion (BIC) 

Another criterion which is often less favorable for models with 
many parameter s is the Bayesian Information Criterio n (BIC) 
dSchwardfl978l) . This criterion, favoured by iLiddld J2004I) . 
weights more heavily additional model parameters than AICc 
for data sets with large N. As AICc, BIC does not require us 
to compare nested models. The BIC consists of computing for 
two models the values: 



= N ln|^| + A,lnA'. 



(6) 



3 On the contrary, the AICc test is less restrictive because it does 
not require nested models. 
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where the parameters have the same meaning as before. As for 
AICc, only the difference of S, is important, not their particular 
value. 

The comparison of the BTF model with the EBTF one gives 
respectively Si = -523.41 and S 2 = -537.34. Applying the 
analogue of Eq. @ for BIC yields that the EBTF model is 1054 
times more likely to be correct than the BTF one. 

3.2.4. Comparison with a quadratic fit 

In order to convince ourself of the strength of the three previous 
comparison methods, we have compared the BTF model with 
a third model, where the data is simply fitted with a quadratic 
function, therefore with the same number of free parameters 
(k = 3) as our model: 

log(M* + M m ) ~ a + b log V c + c log 2 V c . (7) 

Comparing this quadratic model with the BTF model, the F- 
test gives a probability P of only 37.9%. In this case, the weak 
decrease of the rms of the more complex model is not suffi- 
cient to justify the additional parameter c' . The AICc J\ z for the 
two models are respectively -533.79 and -532.51; the BIC S,, 
-523.41 and -518.70 respectively: the BTF model is clearly a 
slightly better model than the quadratic fit. 

Since a quadratic model does not significantly improve the 
fit, we can then conclude that the improvement reached by the 
EBTF model does not result merely from the additional pa- 
rameters but rather from the judicious choice of the fitting law 
Eq. O, directly inspired by physical considerations. 

3.3. Subsamples 

In this Section, we check the robustness of the result obtained 
in Sect. [3]when removing points from the the MSBD data. All 
the subsamples containing all but one or two points are anal- 
ysed. 

For each of the N = 243 subsamples, where in each one a 
single point has been successively removed, we have computed 
the optimal c opt value, as described in Sect. 13. II The histogram 
of those values is displayed at the bottom of Fig. [3] The values 
c opt lie in the interval 2.65 and 3.24 with a mean of 2.98 and a 
dispersion of 0.05. Therefore the optimal c opt is little sensitive 
to the error produced by any single data point. 

The same exercise is repeated by successively removing 
two different points for each distinct combination allowed by 
243 data points. Thus, 243 x 242/2 = 29403 searches for an 
optimal value of c opt have been computed. The found c opt range 
between 2.51 and 3.46 with a mean of 2.98 and a dispersion of 
0.07. 

In the rms of both histogram peaks, we observe that a par- 
ticular removed point contributes to a particularly large de- 
crease of the rms and an increase of c opt . This point corresponds 
to the isolated lower left point of Fig. [5] It is important to em- 
phasize that the conclusions obtained in the previous Section 
are not lessened by this point. On the contrary, removing it 
increases significantly the goodness of the EBTF model with 
respect to the BTF one. 
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Fig. 3. Bottom: Histogram of the optimal c opt values obtained 
by fitting Eq. Q to successively the 243 distinct subsamples 
made of 242 data points. Top: Rms of each subsample fit. 




op(. 

Fig. 4. Bottom: Histogram of the optimal c opt values obtained 
by fitting Eq. Q to each of the distinct 29403 subsamples made 
of 241 data points. Top: Rms of each subsample fit. 

3.4. Independent sample 

In order to check the dependence of the previous results on 
MSBD's particular sample, we have repeated the analysis of 
Sect. l3.2l on another independent sample. We ha ve chosen the 
sample proposed by Sanders & McGaugh (2002) that contains 
76 galaxies, where the luminosities (in either B, R, K or H 
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Fig. 5. Rms of the linear least-squares fit as a function 
of the parameter c (see Eq. for the sub-sample of 
ISanders & McGaughl J2002). The dotted line indicates the rms 
value for c — 1 . 



bands) as well as the HI mass are provided. Among those galax- 
ies, we have removed 29 of them that are contained in MSBD's 
data. They correspond to the UMa sample of Verheiien] J200ll) . 
The average mass-to-light ratio for each band has been cho- 
sen according to MSBD. For the R band, the missing value 
has been interpolated and set to 1.2. The rms of the data fit of 
Eq.^ as a function of c is displayed in Fig. [5] The found opti- 
mal value c opt = 3.42 (rms = 0.256) is slightly higher than for 
MSBD's sample. Beyond this value, the rms grows slowly up 
to c — 16.5 where it is still better than for c — 1 (rms = 0.275). 
The results of the extra sum-of-squares F-Test, AICc and BIC 
tests are displayed in Tabled Those tests compare the BTF fit 
with the EBTF one and with the quadratic fit of Eq. 0. While 
the trend is slightly less pronounced in this new but smaller 
sample, the statistical tests confirm that the EBTF model is bet- 
ter than the BTF one (4 to 9 times). The choice of Eq. (Q) is 
again supported by the fact that the quadratic model (with the 
same number of parameters than the EBTF one) does not im- 
prove the fit. 

3.5. Variable M+/L 

A questionable assumption included in the stellar mass deriva- 
tion in MSBD's sample is that the stellar M+/L ratio is taken 
as constant, while in a sample of spirals it can typically vary by 
a factor of several, more in B than in the NIR bands. 

However, along the spiral sequence the rotational velocity 
is known to correlate with the stellar M*/L ratio: small late 
type spirals typically rotate slower and are bluer (so have a 
lower stellar M* / L ratio) than massive and fast rotating late spi- 
rals. As an illustration, in Broeils' sample (1992), the M*/Lb 



BTF vs EBTF 


BTF vs "quadratic" 


F-Test 


F = 6.81 
P = 1.24% 


F = 2.03 
P= 16% 


AICc 




m\ = -114.72 
3% = -119.08 

Pv/Pia = 8-8 


3{\ = -114.72 
$q = -114.44 
P2.1/P1.2 = 0.87 


BIC 




S, = -109.72 
£ 2 = -112.63 
P 2 ,i/!Pi,2=4.2 


51 = -109.72 

5 2 = -107.99 
P2,i/Pia = 0-42 



Table 1. Statistical comparison of the BTF model with the 
EBTF one, and with the quadratic model, for Sanders & 
McGaugh sub-sample data. 



ratio of 23 bright galaxies spanning the spiral sequence can be 
reasonably represented by a linear regression: 



4.00 log V max - 5.77. 



(8) 



Supposing now that instead of a constant M+/L ratio, one 
would include the above more accurate statistical knowledge in 
the EBTF. The former constant M+/L in MSBD data for each 
galaxy would be replaced by a primed one of the form: 



— - = 1 + d log , 

L \ B 125/ L 



(9) 



where the constant 125 comes from the interval middle value 
of the log of velocities (in km/s) in MSBD's sample; d is an 
additional parameter that might also improve the BTF rms. 
The alternative EBTF fit now reads: 



log 



+ dlog-^jM i , +cM m 



a + blog V c 



(10) 



where here M* is meant not as the true stellar mass, but as 
the quantity provided in MSBD's sample as the stellar mass 
estimate at constant M*/L ratio. 

Repeating the analysis of Sect. 3.1, we search for the 
best pair of parameters (c opt , d ov , t ) that minimizes the rms fit. 
Interestingly, the best fit does indeed require a positive and 
suitable d with regard to the observational data, where 
is correlated with the rotational velocity. This has the effect 
of slightly increasing c. However, the improvement in rms 
is very modest. The best pair is c opt = 3.07, d opt = 0.48, 
for rms = 0.3162, instead of c opt = 2.98, (d opt = 0) and 
rms = 0.3164 for EBTF. The iso-values of rms(c, d) around 
the minimum are strongly elongated, almost straight, parallel 
to the d axis, which already means that the most relevant mini- 
mizing parameter is c. 

With the tests of Sect. 3.2, we can now evaluate quantita- 
tively whether d is a significant additional parameter. For all 
three tests the answer is clearly that d does not bring any sig- 
nificant improvement: for the F-test P - 0.61, for AICc the 
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information ratio is only p2,i/ < Pi,2 = 0.401, and for BIC even 
lower, f 2,1/^1,2 = 0.073. 

The results can be summarized concisely: the principal ef- 
fect of a variable M+/L is to somewhat increase c opt from 2.98 
to 3.07, while suggesting a fairly reasonable M+/L variation 
over the spiral sequences. It especially shows that c is a much 
more relevant parameter than d. 

4. Discussion and conclusion 

An analysis of the baryonic Tully-Fisher relation of MSBD's 
sample shows that a better correlation is obtained by multiply- 
ing the total HI mass by a factor of around 3, suggesting that 
a substantial amount of baryons would remain to be found in 
spirals. 

The optimal factor around 3 is typically lower than the 
caling factor fitting well the total dynamical mass, 5-15 
HoekstraetaDl200lh . There is room for another component 
of a different nature, perhaps non-baryonic and not located 
in the disk, still necessary to explain the rest of the dynam- 
ical mass. This point is supported by the recent work of 
iRevaz & PfennT^erl ll2004 showing that the frequent observa- 
tions of warps in spiral galaxies is the natural result of bending 
instabilities, if, for a galaxy such as the Milky Way, the disk 
mass within 30 kpc represents ~ 60% of the total mass, leaving 
~ 40% for a conventional non-baryonic dark halo. However, 
one should keep in mind that the range of sub-optimal c which 
still improve the original BTF relation rapidly widens toward 
values up to about 11-16, therefore including the full dynam- 
ical mass may still improve the BTF relation. 

In a first step the / L ratio of the stellar component has 
been supposed constant. In a second step it has been allowed to 
vary according to a more realistic trend where M+/L is linear in 
log V c . The statistical tests show that a variable M+/L is com- 
patible with the observations and requires somewhat more dark 
baryons proportional to HI. More clearly the tests demonstrate 
that a variable M*/L ratio is a secondary factor with regard to 
the mass proportional to the HI factor. 

The reason why the Tully-Fisher relation should be a power 
law between the disk baryonic mass and the disk rotational 
velocity remains unclear, because the improvement found on 
the TF relation still does not explain its physical origin; it just 
points toward an explanation based on the disk baryonic con- 
tent instead of either only the stellar or only the detected bary- 
onic content. Dimensional arguments to explain the TF rela- 
tion by internal physical factors in self-regulated gravitating 
disks, instead of initial conditions, were given in Pfenniger 
ill 991 1) where it was pointed out that in pure self-regulated 
disks the gravitational power oc v 5 mt /G that a galaxy can ex- 
change with stars appears comparable to the stellar mechani- 
cal power, about 1/10 of the stellar luminosity in star form- 
ing disks. Remarkably, this estimate of the disk internal power 
provides a fair zero point of the TF relation. The important con- 
sequence is that the mechanical energy output coming mostly 
from massive stars appears sufficient in magnitude to modify 
secularly the global parameters of a disk galaxy. Therefore to 
understand spirals, their internal physics appears more impor- 
tant than the initial conditions of formation. 



This point is consistent with the unified p icture of disk 
galaxies presented in Pfenniaer & Revaz (2004) where bars, 
spirals and warps result from horizontal and vertical instabil- 
ities triggered by the constant energy dissipation. The basic 
properties of these systems may be summarized by an interplay 
of cooling and self-gravitation. In dissipative self-regulated 
gravitating disks, more than the initial conditions, the inter- 
nal physics is crucial to understand their structure and evolu- 
tion. The energy dissipated by radiative cooling, on one hand, 
drives galaxies towards flat disk shapes, but, on the other hand, 
also drives them toward gravitational instability thresholds be- 
yond which horizontal and vertical gravitational instabilities 
are spontaneously triggered, heating the disk both in the hori- 
zontal and vertical directions. This purely dynamical feed-back 
mechanism is supplemented by the heating caused by stellar 
activity jTmmeli et al 12 004). which results from the turbulence 
generated in spiral arms by the larger-scale dynamical heating. 
The galactic disks appear then as tightly self-regulated struc- 
tures that remain in a marginally stable state as long as the 
cooling agent, the gas, exists. This seems to us a suitable start- 
ing point for finding a better physical explanation of the TF 
relation. 
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